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Preface 


Findings presented in this paper are based upon field observations collected as part 
of the Coast and Geodetic Survey’s oceanographic survey program of the central North 
Pacific Ocean. The paper was prepared originally as a thesis in partial fulfillment of the 
requirements for the degree of Master of Science in oceanography at the University of 
Washington, Seattle, Wash., in 1964. As a thesis, the paper was titled ““Deep Circulation 
in a Portion of the North Pacific Ocean During the Summers of 1961, 1962, and 1963.” 
The author, who participated in the North Pacific Ocean survey during 1961, 1962, and 
1963, is grateful for the opportunity afforded him by the Coast and Geodetic Survey 
to pursue graduate study in oceanography at the University of Washington. He would 
like to acknowledge the interest and support shown by the faculty and staff of the Depart- 
ment of Oceanography at the university during preparation of the thesis —particularly 
Dr. Clifford A. Barnes for helpful guidance, Dr. Maurice Rattray, Jr., for critical review, 
Mrs. Monique Rona for machine processing of data, and Mr. Donald R. Doyle for prepara- 
tion of illustrations. Machine processing of data and preparation of drawings were 
supported in part by the U.S. Office of Naval Research (Contract Nonr 477-10, Project 
NR 083012). The paper is published in its original form except for editorial revision of 
format. Illustrations are grouped at the end of the paper. The list of symbols is in 
appendix III. 


iil 


Page 
LEAS EVGA S ss eect Geaereer anc ete HARON o co AEE ATE: BE rie MES TSE et ert arin te ill 
JAI OSTEAG a soooastartamtntl debtany cu Oa SsoRORE CREE Atte esses tN ha eine O aa RL Jb tai OR once em 1 
Pel ntroduction te tenncct cee cc recreate ee Sroreas cea Soret ee Roe ck Come een oe od Cee cent 1 
IPI Objectivessasmery nce ta tae Eee ete eee mee on snetae a tet 1 
NEDAO Dservationsentes noe Cee eae eae es on ec es rea A EEA Ne: 2 
IESE Descriptionvotathekaneaw sos ose te ee echynccrs tac ceeec ono teeeG tee deee macho ioe nenenes 3 
2motate mentaotathes problem sewer eee conan ear ote een eect eee Hoon en neeseoee 5 
MIMO utlinesotsprocedure sce Seacmacscn ese eee eee Ree ee nee 5 
2-2 Assumptions and ap proximationSy.....c2-oseesesse:.tscesececeusceescoceceeenece sees 6 
SMEVE SUL Serena e ener rene AoE Et ae Sant aSCUL move datums the dnia Gatar baths eoatieaee dies 10 
selesblorizomtalllsvelocitiess scence: scence since cues si saciee ec oes ovisice sciajatel oud ersctaolew data does 10 
Sal Velocitie’ smi OOM tes. eadce sects coe oe selaac ae em nak sn tise clgetdde seh aentsadioweeeaty 10 
SRP AVie LOCItIE Su LOO 2a mecrmeiatce ee ose ree conus as oasica ben demise sikiosninae ss we eteee 10 
Dale MME loCItlestIMMLOOS waarcnecaneceee cee mcccerce sate ores chee e accuoneuns aun bee Sarees 1] 
Se MNCL CalavelOCILyeeemeecre tao cre ree eee e en Cece Cat eet eeeiis Ae tate oe eee nna en dese 12 
SeouDepthvornoshorizontallmotionessscccesseste eeeeosasceeeen ec reeeeeeeeciesteneeereeee 12 
3.4 Annual cycle of horizontal velocities..................c.cccececececececencecececeeeacs 12 
At SUMIMaryeanduconclusionssesssencesceseie ner ere eee reece errno eee eee Eee 15 
NGI SIRES soo natcatonu ne coped SOO HER er EG ONC DECOR Seo ECR Ge ER eI= aaNet ein ie 17 
Appendix I. Physical expressions for the problem.....................c.ccecesececeeceeees 19 
Appendixaal EAgnumenicaliexamples-seresrneeseceseeseeercre et cca teecacerc era ostene ceneeaee 27 
AppendixellliPeleistrotesymbolsymecsteseee eee ceee eee ecckecee cece eee eee e eee 35 
iene Smet neous eOG sree scacetits sce errr petites Soe: seiiak ae pM aL rtetaay adhe Pee aero 37-104 
Tables 
TABLE Page 
1. Recent Coast and Geodetic Survey observations in the central North Pacific 
(CYOGE alee ocrccrent soe tis Selec IS IAEA CSCIC a Tene Ee CIs rt PIA fer ae CNP 3 
2. Data printout for section defined by oceanographic stations 63 and 59............. 31 
3. Data printout for section defined by oceanographic stations 29 and 32............. 32 
4. Data printout for section defined by oceanographic stations 63 and 29............... 33 
5. Data printout for section defined by oceanographic stations 59 and 32............... 34 
Illustrations 
FIGURE Page 
IS Locations) of' stationsi and! Se€CtionSis...cn--ecsacee scm sase see ice ce: eisesieiabisisie inte se sletelet iste Be 
2. Stations used in determination of transport...................c eee ee eee eee eee e eee e ees 38 
3. Temperature (°C) Section 1, September 1961..................ececeeeeeneeeeeee ee eeeeees 39 
4. Salinity (°/oo) Section 1, September 1961..................cccceeeeeeeeeeec eee eeeneeeeeeees 40 
5. Dissolved oxygen (ml/L) Section 1, September 1961...................0.e ce ceee sence ees 41 
6. Temperature (°C) Section 2, September 1961............ 0.0... .c ccc cc cence eee eect eee 42 
%. Salinity (/o0) Section 2, September [9G] 2. 2.2. cesses cce cee cne scene ceeesnnaecsnercneees 43 
8. Dissolved oxygen (ml/L) Section 2, September 1961..................0:0ececeeeeeeneees 44 
9. Temperature (°C) Section 3, September 1962...................cccc cece eee eee eee eees 45 
10. Salinity (°/oo) Section 3, September 1962................ 0... cc cece ee ee ee eee eee ee ne nen eees 46 
11. Dissolved oxygen (ml/L) Section 3, September 1962..................0cceseeeeeee eee es 47 


FIGURE 


12: 
13: 
14. 
US: 
16. 
fe 
18. 


19: 


20. 


24. 


20: 


26. 


28. 


29% 


30. 


ols 


32. 


Temperature (°C) Section 4, September 1963...................:ecsssceecneceeeeteceeeees 
Salinity (°/q)psection.4., September 1963.22 sieeceens:---ceecescseraer ees enterctecserss-ee 
Dissolved oxygen (ml/L) Section 4, September 1963.................::cseeeeeeeeeenee ees 
Memperature (¢G)ioection) oy Vay lOO Sreasenentsseeser sense secre teemert ices eee ee 
Salinity ((po) Sections a: May 1963 an ettesc-qunon<scedes ciesahpeeesencsuncces Sehearag tegen 
Dissolved oxyeen (ml/[))Sectiontos May W963 ir esse. scsccsceccseeeteneeenceeeyrsecessse 
Variation with latitude of wind-driven integrated geostrophic transport, Septem- 

ert OG Qos cate ek aoe Sere oe ee see e aeiers aeinc cssleteteinimala nis sieistue ce eolepeecisemeeraresesmeicOemre 
Geopotential anomaly at 4,000 meters referred to 5,000 meters, Sections 2, 3, 


. Topography in dynamic meters, 250-m level adjusted to wind-driven geostrophic 


transport, Ekman divergence considered, September 1961. Drag coefficient 
0.0026; stresses from mean monthly pressureS...................sceeeeeseeeeeeeeeees 
Topography in dynamic meters, 500-m level adjusted to wind-driven geostrophic 
transport, Ekman divergence considered, September 1961. Drag coefficient 
0.0026: stresses from mean monthly pressures................0c eee eee cee ene eee eee es 
Topography in dynamic meters, 1,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, September 1961. Drag 
coefficient 0.0026; stresses from mean monthly pressures............--.:+e1eeee 
Topography in dynamic meters, 2,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, September 1961. Drag 
coefficient 0.0026: stresses from mean monthly pressures.......... +. +.+0eeeee ees 


. Topography in dynamic meters, 4,000-m level adjusted to wind-driven geo- 


strophic transport, Ekman divergence considered, September 1961. Drag 
coefficient 0.0026: stresses from mean monthly pressures...........-+-++eeeeeee es 
Topography in dynamic meters, 5,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, September 1961. Drag 
coefficient 0.0026: stresses from mean monthly pressures..............0++0eeeeees 
Topography in dynamic meters, 5,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence not considered, September 1961. 
Drag coefficient 0.0026; stresses from mean monthly pressures.................- 
Topography in dynamic meters, 2,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence not considered, August 1961. Drag 
coethcient 0.0026; stresses from mean monthly pressures..............:00eeeeeees 
Topography in dynamic meters, 4,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence not considered, August 1961. Drag 
coefhcient 0.0026: stresses from mean monthly pressures..............6.06:ee0e0 
Topography in dynamic meters, 5,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence not considered, August 1961. Drag 
coefficient 0.0026: stresses from mean monthly pressures.............00ee0eee ees 
Topography in dynamic meters, 250-m level adjusted to wind-driven geostrophic 
transport, Ekman divergence considered, September 1962. Drag coefficient 
0.0026: stresses from mean monthly pressures............0.c ce. cee ceeeeeeeeeeneeeens 


vi 


58 


59 


60 


61 


62 


63 


64 


69 


66 


67 


68 


69 


FIGURE 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


42. 


43. 


44. 


45. 


46. 


47. 


48. 


49. 


50. 


Topography in dynamic meters, 500-m level adjusted to wind-driven geostrophic 
transport, Ekman divergence considered, September 1962. Drag coefficient 
0.0026: stresses from mean monthly pressures............... 0.0000 ecceeee ee eens ee ees 

Topography in dynamic meters, 1,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, September 1962. Drag 
coefficient 0.0026; stresses from mean monthly pressures.................000+0+5 

Topography in dynamic meters, 2,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, September 1962. Drag 
coefficient 0.0026: stresses from mean monthly pressures................0.0.500+ 

Topography in dynamic meters, 4,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, September 1962. Drag 
coefficient 0.0026: stresses from mean monthly pressures....................605- 

Topography in dynamic meters, 5,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, September 1962. Drag 
coefficient 0.0026; stresses from mean monthly pressures....................-5-+ 

Topography in dynamic meters, 5,000-m level adjusted to wind-driven geostro- 
phic transport, Ekman divergence not considered, September 1962. Drag 
coefficient 0.0026: meaned stresses from twice-daily reports of pressure...... 

Topography in dynamic meters, 2,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence not considered, September 1962. Drag 
coefficient 0.0012: meaned stresses from twice-daily reports of pressure...... 

Topography in dynamic meters, 5,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence not considered, September 1962. Drag 


Topography in dynamic meters, 250-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, May 1963. Drag coefh- 
cient 0.0026; stresses from mean monthly pressures..............0:0seeee serene nes 

Topography in dynamic meters, 500-m level adjusted to wind-driven geostrophic 
transport, Ekman divergence considered, May 1963. Drag coefficient 0.0026; 
stresses from mean monthly pressures..............2:seeceeeeeceec escent eee ee eee eeees 

Topography in dynamic meters, 1,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, May 1963. Drag coefh- 
cient 0.0026; stresses from mean monthly pressures..............0.0e0eeeee eee ees 

Topography in dynamic meters, 2,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, May 1963. Drag co- 
efficient 0.0026; stresses from mean monthly pressures..................0.0202000 

Topography in dynamic meters, 4,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, May 1963. Drag co- 
efficient 0.0026; stresses from mean monthly pressures..................e0eeee ee 

Topography in dynamic meters, 5,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence considered, May 1963. Drag coefh- 
cient 0.0026; stresses from mean monthly pressures..................0.00eeeee eee 

Topography in dynamic meters, 5,000-m level adjusted to wind-driven geo- 
strophic transport, Ekman divergence not considered, May 1963. Drag 
coefficient 0.0026; stresses from mean monthly pressures....................0006 

Topography in dynamic meters, 250-m level adjusted to wind-driven geostrophic 
transport, Ekman divergence considered, May 1963. Drag coefficient 
0.0012: meaned stresses from twice-daily reports of pressure...............60005- 

Topography in dynamic meters, 500-m le ‘el adjusted to wind-driven geostrophic 
transport, Ekman divergence considered, May 1963. Drag coefficient 
0.0012; meaned stresses from twice-daily reports of pressure...................5- 


vil 


Page 


70 


71 


1 


74 


75 


76 


il 


78 


79 


80 


81 


82 


83 


84 


85 


86 


FIGURE Page 
51. Topography in dynamic meters, 1,000-m level adjusted to wind-driven geo- 


Spy 


34. 


strophic transport, Ekman divergence considered, May 1963. Drag coefh- 

cient 0.0012: meaned stresses from twice-daily reports of pressure............. 87 
Topography in dynamic meters, 2,000-m level adjusted to wind-driven geo- 

strophic transport, Ekman divergence considered, May 1963. Drag coefh- 

cient 0.0012: meaned stresses from twice-daily reports of pressure............. 88 
Topography in dynamic meters, 4,000-m level adjusted to wind-driven geo- 

strophic transport, Ekman divergence considered, May 1963. Drag coefh- 

cient 0.0012; meaned stresses from twice-daily reports of pressure............. 89 
Topography in dynamic meters, 5,000-m level adjusted to wind-driven geo- 

strophic transport, Ekman divergence considered, May 1963. Drag coefh- 


cient 0.0012: meaned stresses from twice-daily reports of pressure............. 90 

. Variation of vertical velocity with depth at indicated latitudes along approxi- 
mately 1 70> Wiese cats-sacas eietias san eatedtacis ceca seman schtrocnmatemeu sacs anaeeayeeneterts 9] 
. Depthtotnorzonalamotion lO0m Wg LOO. eum seeeiense cient sreeereemeeetter engines 92 
2 Denth oto zonal motions 1 (ia W 096i eseeneent cen emteaaan ce seer rere reemersrtrsererre 93 
_ Depth of no meridional/motion, 169° W, 196. 222. .c cee cas sac ec es seenccesteeaecees 94. 
SDepthvot novzonalemo tions lO 0s eWeOO2 meres ccet ceverterseee sae Ceca ae smee secs et eer: 95 
. Depth of no zonal motion, 1657 W. May 1963.00 3.csc. cscs. 2. seemencnnncts scree nestereen 96 
» Depihyot no: zonalimotions 180 2W saya 963 secrete yee 2.2: - eases esmencemereeeeene 97 
2. Depth of no mendional’ motion, 1727 W, May 1963...0... 2.s.cccc aes -2c ends eeeeeeeenes 98 
» Annual variations in zonal velocity at 4000 ines. 252s. essen seer se ses sec ce steels 99 
4. Variation with latitude of integrated total transport, September 1961............ 100 
. Variation with latitude of zonal Ekman transport, September 1961............... 101 
. Variation with latitude of zonal Ekman transport, September 1961............... 102 
. Variation with longitude of meridional Ekman transport, September 1961......... 103 
. Variation with latitude of integrated geostrophic transport, September 1961... 104 


vill 


Deep Circulation, Central North Pacific Ocean: 
1961, 1962, 1963 


Abstract. Methods are reviewed for obtaining velocity reference values such that the absolute flow in the 
deep sea can be described in terms of these reference values and the dynamic height computations. The 
method of Murty and Rattray (1962) for obtaining absolute flow by adjusting the transports from dynamic height 
computations to analogous values of transport determined from the curl of the wind stress was examined and 
used for calculating velocities for a portion of the central North Pacific Ocean between the Hawaiian and Aleutian 
Islands from observations made during the summers of 1961, 1962, and 1963. 

A predominantly eastward flow was found at depths of 250 and 500 meters. This flow is similar to the 
accepted pattern of flow in this area. Over much of the area velocities at 4,000 and 5,000 meters were on the 
order of 1 cm/sec and in a direction opposite to those near the surface. The flow near the bottom tends to con- 
form to salient features of the bathymetry. By comparison with flow inferred from the distribution of tem- 
perature and dissolved oxygen, it is concluded that the flow at depths of 4,000 meters and more computed for 
the summer months does not reflect the net annual circulation, but is a periodic barotropic flow. 

Errors in the assumptions in the formulation of the equations used for computing the flow, in the inadequacies 
in tabulations of wind-driven transport, and in observational errors are examined. 


1. Introduction 
1.1 Objectives 


The method of computing relative velocities by the dynamic heights method has be- 
come standard in oceanography. Dynamic height computations can give a realistic 
representation of the absolute velocity field in respect to the earth if proper allowances 
can be made for wind drift and if the relative velocities can be referred to absolute veloci- 
ties at some level. However, the estimation of absolute velocities for reference has often 
lacked rigor. 

The main objective of this paper is to describe as realistically as possible, the 
circulation in a portion of the North Pacific Ocean. A method is developed which relates 
transports calculated from oceanographic station data to geostrophic transports deter- 
mined from the wind-stress field to find values of the velocity near the bottom. These 
velocities near the bottom are then used as references for the dynamic height velocities. 
The method is applied using three recent series of oceanographic observations and tabu- 
lations of wind-driven geostrophic transports from Fofonoff (1960). 

Defant (1961) reviewed methods of defining and delineating a reference surface for 
dynamic topography. He cited the method of determining absolute from relative values 
by comparing direct current measurements with the relative velocities. However, he 
pointed out that the difficulty in making a sufficient number of direct measurements, each 
long enough to define the mean field of motion, drastically limits the use of this method. 

In practice, a reference surface is selected arbitrarily near the ocean bottom or as 
deep as the observations allow. Surface and near-surface flow generally can be re- 
alistically estimated because velocities at the reference depth are usually negligible in 
comparison to velocities near the surface. However, this arbitrary selection of a reference 
level can result in relatively large errors in computing velocities of the deeper currents, 


and in computations of transport. Generally, the depth of zero velocity will not be a level 
surface, but will vary from place to place. In some early attempts to define this variable- 
depth reference surface, it was placed at the depth of minimum dissolved oxygen. How- 
ever, it is difficult to accept the assumption that the minimum of dissolved oxygen is a 
consequence of no horizontal advection. Further, velocities obtained using this reference 
surface were unrealistic. 


Defant calculated the depth of the reference surface by a method which is implicit 
in the dynamic height observations. He used his method to construct a chart of the 
depth of the reference surface in the Atlantic Ocean. A plot was made of pressure (or 
depth) versus the difference in dynamic heights for each pair of adjacent stations. The 
plotted curve represents the velocity component normal to the section. Any part of the 
plotted curve that is parallel to the pressure coordinate represents a depth interval of 
constant velocity. Defant reasoned that the velocity is more likely to be zero in this layer 
of no shear stress than at any other level. In constructing his chart of the depth of this 
reference surface, he made the additional assumption that the horizontal flow parallel 
to the dynamic height section is also zero at this level. Although this method of selecting 
a reference surface is subjective, its use has resulted in consistently realistic velocity 
fields. This is essentially the method used by Russian investigators working with 1958— 
1959 Vityaz data from the North Pacific Ocean (Chekotillo, 1961). 


Stommel (1956) described a method for determining the depth of no meridional motion 
by equating the divergence of the Ekman transport with the divergence of the geostrophic 
transport. Using the concept introduced by Fofonoff (1962) that the total transport can 
be separated into a barotropic mode, a baroclinic mode, and an Ekman transport, Stom- 
mel’s method can be expressed as a relationship between these three modes: 


Lge 
PmVUb = i (V _ Vr ee Vz) 


where the V’s represent mass transport per unit width, wm is the velocity at the bottom, 
h is the total depth, pm is the mean density of the water column, and the subscripts » and g 
denote Ekman and baroclinic modes, respectively. Stommel did not consider the effects 
of bottom topography on the divergence of the barotropic mode of transport. 

Murty and Rattray (1962), Fofonoff (1961), Favorite (1961), and Dodimead, Favorite, 
and Hirano (1962) have compared baroclinic transports with wind-driven transports. 
The method used by Murty and Rattray, in which the baroclinic transport is compared 
with the wind-driven geostrophic transport to obtain a velocity near the bottom, is used 
herein. The method is developed in appendix I and a numerical example is given in ap- 
pendix II. 

Russian workers have determined deep currents by similar methods. Koshlyakov 
(1961) reported computations in the North Pacific Ocean based on arguments similar to 
Stommel’s, but included the effects of variations in depth by means of an expression 
similar to: 


7 = OTsy OTs 

Re i (7 is) 
PmVb a(fh) 
dy 


In this expression the new terms Ts and 7s, are wind stresses at the surface, f is the rate 
of change of the Coriolis parameter with latitude, and f is the Coriolis parameter. This 
expression is valid in areas where the product of the zonal component of velocity at the 
bottom and the zonal component of the bottom slope is much smaller than the product 
of the meridional components of the same variables. 
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1.2 Observations 


In 1961, the Coast and Geodetic Survey of the United States Department of Commerce 
began field work on a survey of that part of the North Pacific Ocean between the Hawaiian 
Islands on the south, the Aleutian Islands on the north, and (approximately) the 158° W 
and 180° meridians on the east and west, respectively. Oceanographic stations in this 
program included observations of temperature by reversing thermometers, salinity by 
salinometer, and dissolved oxygen by Winkler titration (table 1 and fig. 1). Only those 
stations deeper than 2,000 meters (fig. 2) were useful for the present study. The field 
observations are summarized graphically in figures 3 through 17. 


TABLE 1.—Recent Coast and Geodetic Survey observations in the central North Pacific 


Ocean 
Number | National Oceanographic 
Year Dates Vessel Sections of deep Data Center reference 
| stations number 

1961 10 Sept.-12 Oct. Pioneer 160 °W. long. 6 31099 
177 °W. long. 6 31099 
1962 22 Aug.—25 Sept. Pioneer 160 °W. long. 10 31099 
1963 6-29 May Pioneer 165 °W. long. LOR WA ep. eaten: 
6-29 May Surveyor 180 ° long. Cmerlly E ae te 


Values for the variables at standard depths and computed values of sigma-t, the 
specific volume anomaly, the geopotential anomaly, and the potential energy anomaly 
(Fofonoff, 1962), were determined by means of an interpolation and computation program 
prepared by the Applied Mathematics staff of the Department of Oceanography, University 
of Washington. The computations were made on the University’s IBM 709 computer. 
The 1961 and 1962 observations were also submitted to the National Oceanographic Data 
Center, where values for standard depths were again calculated. The values reported 
by the National Oceanographic Data Center disagree slightly with those from the Uni- 
versity because of minor variations in the two interpolation and computation programs. 


1.3 Description of the Area 


The observations listed above are of considerable interest in supporting recent de- 
criptions of the fields of oceanographic variables in the central North Pacific Ocean. 
Since 1955, co-operative efforts such as those summarized by the NORPAC Committee 
(1960a, b), observations resulting from International Geophysical Year cruises, programs 
of the International North Pacific Fisheries Commission and of the Bureau of Commercial 
Fisheries in connection with its Pacific Oceanic Fisheries Investigations, have provided 
a basis for detailed descriptions of the area by previous investigators. The present 
investigation confirms these descriptions of temperature, salinity and oxygen distributions 
(figs. 3 to 17 inclusive), and amplifies certain details. However, the primary purpose of 
this investigation is to describe the circulation and characteristics of the deep water. 

Most oceanographic descriptions of this area treat the upper 1,000 meters because 
that stratum includes the more striking parts of the structure, and relatively few data 
are available from greater depths. Hence, most descriptions of the circulation deal with 
surface velocities or are limited to velocities and transports above 1,000 meters. Reid 
(1961) described the surface geostrophic flow for the Pacific Ocean. His calculations 
were based on a level of no net motion at 1,000 decibars. The NORPAC Atlas (1960a) 
covers in greater detail the geostrophic flow in the North Pacific Ocean. Permanent 
features of the surface flow in the central north Pacific Ocean are: the Alaskan Stream and 


3 


Alaskan Stream Extension, which form a relatively intense westerly flow south of and 
closely confined to the Aleutian Island Arc. A zone of easterly flow lies between 30° N 
and 50° N. Dodimead, Favorite, and Hirano (in press) designate parts of this easterly 
flow as the Subarctic Current, made up in large part of Oyashio waters: the Westwind 
Drift, originating where parts of the Oyashio and Kuroshio meet and turn eastward: and 
the North Pacific Current, formed by the main flow of the Kuroshio which turns easterly 
at about 35° N. From 30° N to the southern limits of the present study area the flow is 
weak, as would be expected near the center of the anticyclonic flow in the subtropical 
north Pacific Ocean. A westerly flow persists at the southern limit (23°30'N) of the 
present investigation. 

Strong gradients of temperature and salinity are found across the Westwind Drift. 
This band of strong gradients, extending nearly across the north Pacific Ocean is called 
the Polar Front. Fleming (1955, 1958), Uda (1963), and Dodimead, Favorite and Hirano 
(in press) divide the North Pacific Ocean into natural regions based on the current systems 
and the associated characteristics of temperature, salinity and dissolved oxygen. The 
southern edge of the Polar Front is the boundary between the Subarctic Region in the north 
and the Subtropic Region in the south. Dodimead and coworkers define the boundary 
between the two regions as the zone where the 34°/,, isohaline is nearly vertical from the 
surface to a depth of 300 to 500 meters (figs. 4, 7, 10, 13, and 16). In the present study area, 
the boundary occurs at about 40° N. 

Above 1,000 meters, the density structure in the Subarctic Region is dominated by 
salinity effects. Seasonal variations in the salinity structure are exhibited only in the zone 
extending from the surface to about 100 meters. Below this zone a permanent halocline, 
in which the salinity increases rapidly, occurs in the 100- to 200-meter depth interval. 
Below the halocline, the salinity gradually increases with depth. Fleming (1958), and 
Tully and Barber (1960) have described the salinity structure in detail. The negligible 
seasonal variation in all but the upper 100 meters of the water column is explicit in their 
descriptions. The temperature and salinity structures in the Subarctic Region are subject 
to pronounced seasonal variations through the same depth interval. A temperature mini- 
mum is often present in the depth interval of the halocline as shown in the northern end 
of the sections (figs. 6,9, and 12). The salinity distribution dominates the density struc- 
ture and maintains stability despite the temperature inversion. 

The halocline becomes progressively less pronounced southward through the Polar 
Front, and finally disappears. 

In the Subtropic Region there is no permanent halocline. Here the density structure 
in the upper 1,000 meters is dominated by a strong thermocline extending to depths of 600 
to 800 meters (figs. 3, 6, 9, 12, and 15). A salinity minimum occurs at depths of 500 to 
600 meters (figs. 4, 7, 10, 13, and 16). Seasonal effects in the surface layers in the Sub- 
tropic Region are less pronounced than in the Subarctic Region. 

The division of the North Pacific Ocean into regions on the basis of the temperature 
and salinity distributions applies between the surface and 1,000 meters. Below 1,000 
meters, temperature and salinity structures are similar throughout the area. 

The dissolved oxygen structure is characterized by a minimum occuring between 600 
and 1,200 meters. It occurs at about 600 meters immediately south of the Aleutian Island 
Arc, descends to about 1,200 meters in the region of the Polar Front, and ascends through 
the Subtropic Region to about 800 meters near 23°30'N. <A noteworthy feature of the 
meridional sections of oxygen concentration (figs. 5, 8, 11, 14, and 17) is the thickening 
of the zone of low oxygen values in the northern part of the area. From 46°N to 52°N, 
the thickness of the oxygen zone, of concentrations less than 1 milliliter/liter, is twice 
that south of 40°N. 

Figures 19 through 54 show the results of the dynamic computations for the sections 
investigated in this study. These will be considered in the text which follows. 
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2. Statement of the Problem 


2.1 Outline of Procedure 


The basic premise is that the mass transport calculations tabulated by Fofonoff (1960) 
describe the total transport, and that they can be used as reference values for transport 
calculations by the dynamic heights method. The procedure for determining these 
reference values is as follows: 

a. The geostrophic mode of the total transport is extracted from Fofonoff’s tabula- 
tion for the section separating two oceanographic stations. 

b. Transport through the same section is computed from the dynamic heights at 
the two stations. These dynamic height calculations are based on a level at 
4,000 or 5,000 meters, whichever more nearly equals the mean depth of the 
section. 

c. The transport determined in b (as described by the dynamic heights) is the 
total baroclinic transport (provided the flow due to internal structure vanishes 
at the bottom) and the difference between the transports determined by methods 
a and b must be the barotropic transport. This barotropic transport divided 
by the water density and the area of the section gives the barotropic velocity. 

d. Combining the barotropic velocity with baroclinic velocities determined from 
the dynamic height calculations results in a description of the total geostrophic 
flow through the section. If the assumption in c (that baroclinic flow is zero at 
the bottom) is not valid, then the division of the flow into barotropic and baro- 
clinic modes is indeterminate; but the description of the total geostrophic flow 
and its variation with depth is still valid. 


To establish a basis in theory for the method, to determine the limitations of (and the 
assumptions made in) Fofonoff’s (1960) transport calculations, and to justify using his re- 
sults in the present manner, it is necessary to develop applicable expressions for transport. 
This has been done in appendix I. Based on Fofonoff’s (1962) treatment of Munk’s (1950) 
wind-driven flow formulation, the expressions developed for barotropic velocities u, and 
Vp are: 


bb = 
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(35a, b) 


In (35a), y is a transport function for total transport, and so the first term in the numerator 
on the right expresses the total transport through a meridional section. Zonal Ekman 
transport per unit length along the meridian is denoted Ug, R is the radius of the earth, 
and Ad is a change in latitude, so the second term is Ekman transport through the same 
section. The final term is baroclinic transport, expressed as the difference in potential 
energy anomalies, x, across the section, divided by the Coriolis parameter, f The de- 
nominator on the right side is the mean density, Pm, times the area of the section, hARA@. 
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2.2 Assumptions and Approximations 


The validity of results from the approach used here is dependent on the assumptions 
(a) that the wind-driven geostrophic transport tabulated by Fofonoff (1960) is accurate 
in the region of interest: (b) that acceleration terms, bottom friction terms and horizontal 
friction terms can be neglected in developing the expressions for wind-driven transport 
(p. 21, appendix I); and (c) that convergence due to the effect of bottom topography can 
be neglected. The results will also depend on (d) the accuracy, precision, number and 
distribution of the oceanographic observations. 

(a) Fofonoff (1960, p. 3), in considering the applicability of his computations of wind 
stress at the surface, stated, “It should be remembered that the transports are very sen- 
sitive to the proportionality factors used in relating geostrophic to surface wind and 
surface wind to wind stress. Numerical equivalence of computed and observed transports 
is not anticipated.” His calculations of surface stress are based on the equation 


Ce tol PaC pw’, 
where 7; is the surface stress, Pa is the air density, Cp is the drag coefficient, and w is the 
wind speed. It is assumed that the wind speed and stress magnitude are functions of 
the pressure difference, Ap, between positions so that 


Ts © w? & (Ap). 


Additional error is introduced in Fofonoff’s computed stresses by using the square of 


monthly means of pressure differences, (Ap)’, rather than computing the mean stress from 


the actual pressure (that is, satisfying the relationship 7; * Ap’). The error arising from 
calculation of stresses from a monthly mean of pressure differences rather than using 
the monthly mean of stresses estimated from continuous reports of pressure is propor- 


tional to the variance of the pressure differences. Since Ap? is larger than (Ap), the less 
rigorous relationship leads to stresses that are too small. 

Fofonoff (1960) used a value of 0.0026 for the drag coefficient in his tabulation of 
wind-driven transport. Deacon and Webb (1962) have summarized recent determinations 
and find a range of values from about 0.0010 to 0.0015. Because of the linear relation- 
ships between the drag coefhcient and Ekman and total transport, this discrepancy in- 
creases Fofonoff’s values by a factor of about 2. 


Graphic comparisons have been made to assess the effects of these two discrepancies. 
Figure 18 shows the alternate determinations for the variation with latitude of integrated 
geostrophic transport in September 1962 (Section 3). The three determinations are 
based on: 

~ (i) stresses obtained by using a drag coefficient of 0.0026 and monthly means of 

pressure differences; 
(ii) a drag coefficient of 0.0026 and a monthly mean of values estimated from twice- 
daily pressure reports; and, 
(iii) a drag coefficient of 0.0012 and a monthly mean of the values estimated from 
twice-daily pressure reports. 
The use of the higher drag coefficient (of 0.0026) rather than the more nearly correct 
value of 0.0012 leads to high transport values; but the use of monthly mean pressures 
(rather than means from twice-daily values) gives low transports. Accordingly, errors in 
method (i) are partly compensated and the results obtained agreed well with method (iii); 
whereas method (ii) gives erroneously high results in all cases. Figures 38, 39 and 41 
show for the 5,000 meter level the dynamic topography resulting from the September 1962, 
observations adjusted to these three techniques of wind-driven transport computation. 
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The effects of adjusting the dynamic topography to the geostrophic transport by 
methods (i) and (iii) can be seen by comparing figures 42 to 47 (drag coefficient of 0.0026 
and monthly mean pressures) with figures 49 to 54 (drag coefficient 0.0012 and twice-daily 
mean pressures). These figures show the same general pattern and similar slopes with 
discrepancies of only | to 2 dynamic centimeters in a few locations. 

(6) The calculations for wind-driven transport for each month are based on the stress- 
distribution for that month’s mean pressure distribution, without consideration of initial 
conditions or the inertial response characteristics of the ocean. For the assumption to 
hold that acceleration terms are negligible, the ocean system must be essentially in equi- 
librium with wind-imposed forces during each month considered. Intuitively, a response 
time shorter than one month would seem realistic for the barotropic mode, but intuition 
as well as analysis of observations precludes such short response times for the baroclinic 
mode. Veronis and Stommel (1956), in their treatment of variable wind stress on a two- 
layered ocean without lateral boundaries, showed that for periods of the order of one month, 
acceleration terms were unimportant in both the barotropic and baroclinic modes. An 
explanation for this is that the barotropic mode has a response time much shorter than one 
month and so it is at all times essentially in equilibrium with the driving force, while the 
baroclinic mode primarily responds to forces with periods of the order of one year or more. 
Presumably fluctuations of wind stress for periods shorter than one month will give in- 
significant accelerations. 

In addition to wind stress, bottom and lateral friction must also be considered. Proud- 
man (1953) calculated bottom stress using the square law 


= 2 
Th — ChPbu, 


in which the drag coefficient, c;, is taken as 0.0025, p,, the density of water at the bottom 
is roughly 1 gm/cm?*, and uy, the velocity in cm/sec, is measured one meter above the bot- 
tom. In the present investigation, deep velocities on the order of 1 cm/sec were found. 
Under these conditions, bottom stresses would be of the order of 0.002 dynes/em? and 
negligible compared to surface stresses. 

Munk (1950) developed expressions for wind-driven transport in a central region of 
the ocean in which horizontal friction terms were neglected. These expressions result 
in a distribution of transport conforming to the distribution of properties, and thus support 
the assumption that horizontal friction is negligible in central oceanic regions. 

(c) The assumption that divergence due to the effect of bottom topography is neg- 
ligible in the area of study is examined here. Sverdrup (1947), Munk (1950), and Fofonoff 
(1962) all specified that the transport was made up of only baroclinic and Ekman modes; 
therefore, the treatment of the bottom convergence problem was not pertinent. If equa- 
tions (22), (23), and (24) from appendix I are added and the term for bottom effects retained: 


OT sy OTs P F d0Zy OZp 
BV gt+VetVi)= oe le ea 


The term for bottom effects can be expressed as 
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) 
where @ is the angle between the bottom velocity and gradient vectors. Substituting and 
rearranging, 
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The meridional barotropic transport can also be approximately expressed as 
V, = hvppm. 


As noted in the discussion of the accuracy and precision of observations in section (d) 
which follows, velocities are significant to within about 0.6 cm/sec. Representative 
values of the depth, h, and the mean density, pm, are 500,000 cm and 1 gm/cm?, respec- 
tively, so the resulting difference in barotropic transport for a velocity difference of 0.6 
cm/sec would be 


V,~ 3-10° gm cm™! sec}, 


and terms on the right side of the expression for barotropic transport of order 104 gm em7! 
sec! and less could be neglected. Total transport and barotropic transport terms are of 
the order 10 gm cm7! sec7!. The Ekman transport is of the order 104 or 10° gm cm! 
SeGua: 

In examining the data from September 1961 (fig. 28), it is apparent that the most critical 
combination of bottom velocity and bottom slope occurs at the Hawaiian Island Rise. 
Between stations 25 and 26, the bottom velocity is about 2 cm/sec. The bottom slope of 
2 - 10-* was computed from the smoothed isobaths on H.O. Chart 5486 and is in general 
agreement with that computed from station soundings. The angle 6 which the stream 
lines make with a line normal to the isobaths is uncertain because of the lack of detail 
of bottom topography and the general character of the velocity determinations. It is ap- 
parent, however, that the bottom water tends to flow along the bottom contours, a tendency 
particularly evident in the northern and southwestern portions of the area observed (figs 
27 and 28). Assuming a value of 6 of 85° (suggested in the present data by the close ap- 
proach of the direction of the bottom currents to the isobaths), bottom velocity as 2 cm/sec, 


f as 3X 108 for 25° N, and p=1 gm/cm®*, the magnitude of the final term in the expres- 
sion for barotropic transport is approximately 10° gm cm~! sec-!. This magnitude in a 
depth of 5,000 meters results in a constant velocity of 0.2 cm/sec —one third the significant 
value of 0.6 cm/sec (see below). 

(d) Since the ratios of the baroclinic to barotropic modes within the flow are computed 
from oceanographic observations, the accuracy and the precision of these observations 
must be examined for their influence on the velocity at any level. Velocity at the bottom 
is of primary interest. If values of +0.01°C and +0.01°/o9 are accepted as the precision 
of temperature and salinity measurements and the difference between stations is such as to 
give the maximum error in specific volume anomaly for the interval of 0 to 1,000 decibars, 
an error of £0.01 dynamic meter will arise. Wooster and Taft (1958), using a statistical 
approach to this problem, found that the standard error of the difference in geopotential 
anomaly (0 over 1,000 db) between two stations is 0.0056 dynamic meter and conclude that 
the measurement error of this difference is +0.011 dynamic meter (26). For a pair of 
stations at 30° latitude separated by 250 km, the error in the computed geostrophic velocity 
is +0.6 cm/sec. These values of anomaly and velocity must be kept in mind when asses- 
sing the significance of the horizontal velocities determined in this program. The method 
used, equating baroclinic transport plus barotropic transport to wind-driven geostrophic 
transport, partially compensates for errors in the observations. Any observational error 
leading to a baroclinic transport too large will result in barotropic transport too small by 
the same amount. Thus the barotropic velocity, constant with depth, will be slightly too 
small. The error in total computed velocity at the depth of the observational error will be 
slightly less than that in the baroclinic velocity and almost negligible at other depths. 
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The significant value of geopotential anomaly (0.01 dynamic meter) must also be con- 
sidered in assessing the significance of computed vertical velocities. Since the flow 
determined here is geostrophic, the variation of vertical velocity, w, with depth can be 
approached through the expression 


ow Bp 
a 
4 if 

where vis the meridional component of total velocity and 6 and fare as previously defined. 
In the present study area the horizontal flow is mainly zonal; meridional velocities are 
usually of the same order as the level of significance determined for horizontal velocities. 


dw : : ; 
So the values of ae obtained by the above expression are questionable. 


To help establish the consistency of the observations, and also to determine the validity 
of the proposition that baroclinic velocity vanishes at the bottom, plots were made of the 
geopotential anomaly as a function of latitude. Figures 19 and 20 show the geopotential 
anomaly at 4,000 meters referred to 5,000 meters; figures 21 and 22 are of the anomaly 
at 3,000 meters referred to 5,000 meters. At 4,000 meters, the greatest range of the 
anomaly is 0.02 dynamic meter, only twice the measurement error. This does not permit 
the assumption that baroclinic flow in the interval from 5,000 to 4,000 meters can be 
ignored, but it does indicate the relative unimportance of this deep internal structure. 
The geopotential anomaly at 3,000 meters referred to 5,000 meters varies more with time 
and place than does the anomaly at 4,000 meters. Effects of water structure in the 3,000- 
to 4,000-meter interval must certainly be included in the determinations. 
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3. Results 


3.1 Horizontal Velocities 


Dynamic heights were calculated and adjusted to give transports equal to those de- 
termined from the curl of the wind stress (see appendix IJ for anumerical example). These 
adjusted heights are shown as lines of equal geopotential through the area of the investi- 
gation (figs. 23 through 54). The selected contour interval is 0.01 dynamic meter, a value 
which is uncertain in the upper levels where the paucity of stations must result in excessive 
smoothing of the dynamic topography. But this interval is used throughout to facilitate 
comparison of all levels. 


3.1.1 Velocities in 1961.—In figures 23 through 28, the wind-driven geostrophic trans- 
port for September 1961 was used to obtain reference values for dynamic height calcula- 
tions based on oceanographic observations made during that month. Geopotential 
anomalies at the 250, 500, 1,000, 2,000, 4,000, and 5,000 meter levels are presented. The 
geostrophic transport for these figures was determined from the difference between wind- 
driven total transport and Ekman transport. Figure 29 shows the geopotential anomalies 
at 5,000 meters calculated from the September 1961 oceanographic data and the wind- 
driven geostrophic transport which was computed from the tabulated geostrophic transport 
function. Figures 30 through 32 show the results of a comparison of the flow computed 
from the geostrophic transport function in August 1961. 

The flow pattern was generally zonal, especially in the upper layers. North of about 
40° N, the direction of flow was easterly at 250, 500, and 1,000 meters, but reversed to 
westerly at 2,000, 4,000, and 5,000 meters. South of 30° N, a significant southerly flow 
occurred at 250 and 500 meters. At 2,000 meters, and below, the flow conformed to the 
trend of the bottom contours along the Hawaiian Island Rise. The deep velocities were 
westerly north of the rise and easterly south of it. Figures 28 and 29 can be used to com- 
pare the consequences of using the wind-driven flow with and without consideration of 
Ekman divergence to obtain reference values for dynamic heights. No significant differ- 
ences are noted. An analysis of the oceanographic observations used in these determina- 
tions (or of other observations in the same area) shows the water structure leading to 
baroclinic transport to vary only slightly. A comparison between September 1961 oceano- 
graphic data adjusted for transports from the wind stresses of September (figs. 26, 27, and 
28) and then of August (figs. 30, 31, and 32) gives an indication of the variations in the 
barotropic velocity. 

The flow at 4,000 and at 5,000 meters tended to conform to the isobaths. This ten- 
dency was evident both at the Hawaiian Island Rise and near the Aleutian Trench. 


3.1.2 Velocities in 1962.—The oceanographic observations in 1962 comprise a single 
section along the 160°W meridian. By adjusting the values of the dynamic heights to 
yield transports equal to the geostrophic transports arising from different wind stress dis- 
tributions, an evaluation is made of the effects of (a) estimating wind stress from a monthly: 
mean of atmospheric pressure distributions contrasted to estimating it from twice-daily 
reports of pressure distribution and taking a monthly mean of these values, and (b) using 
different drag coefficients in combination with the two alternate methods of finding the 
mean monthly stress values. Figures 33 through 38 show the geopotential topography 
adjusted to geostrophic transport arising from wind stresses based on a drag coefficient 
of 0.0026 and the mean pressure distribution for September 1961. Figure 39 is based on 
the same drag coefficient but the average of two daily reports of atmospheric pressure 
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were used in computing stresses. Figures 40 and 41 are based on a drag coefficient of 
0.0012 and twice-daily reports of atmospheric pressure. 


The flow patterns as computed vary with the method of calculating stresses and 
particularly marked differences in these computed patterns occur in high latitudes. Com- 
parison of figure 38 with figure 39 shows that alternative methods of computing mean 
stresses give apparent velocities in opposite directions north of 45°N. However, when 
both the drag coefficient and the method of obtaining mean stresses were changed (figs. 
40 and 41), the resulting velocities at equal levels were but little different up to 45°N from 
those indicated in figures 36 and 38. It should be noted that Ekman divergence is not 
considered in figures 36 and 38, and its significance is discussed in a following section. 
The use of 0.0026 as the drag coefficient resulted in values for stress that apparently were 
too large, but calculating stress values from a monthly mean of pressure distribution rather 
than estimating stresses from twice-daily reports of pressure and then taking a monthly 
mean led to stress values that apparently were too small, and these discrepancies tend to 


cancel one another. 
3.1.3 Velocities in 1963.—The oceanographic data for 1963 are presented in the same 


manner as the 1961 data. Figures 42 through 47 are of geopotential topography adjusted 
to wind-driven geostrophic transport with consideration of Ekman divergence. The 
wind-driven transport is based on a drag coefficient of 0.0026, and the mean pressure 
distribution for May 1963 is used. As in the 1961 results, flow was mainly zonal. North 
of 40° N, the flow at 250, 500, and 1,000 meters was easterly. Velocities decreased from 
about 3 cm/sec at 250 meters to less than 1 cm/sec at 1,000 meters. From 40° to 30° N, 
flow was southerly. At 500 meters, flow was easterly or northeasterly north of 30° N: it 
was southwesterly south of 30° N. At the 1,000-meter level, the flow south of 40° N was 
poorly defined in the 1963 data, though there were indications of a cyclonic flow centered 
at about 35° N, 180°. South of 28° N, there was a weak westerly set. At 2,000 meters, 
flow was weak over most of the area and there is some indication of an anticyclonic cir- 
culation centered near 35° N, 180°. The patterns of the flow at 4,000 and 5,000 meters 
were similar to each other. North of 45° N, velocities at these depths were of the order 
of 1 cm/sec westerly, and were the largest found at or below 4,000 meters in 1963. Be- 
tween 45° and 40° N, the flow was easterly, shifting to southerly south of 40° N. 

Flow at the deep levels tended to conform to the salient features of the bottom topog- 
raphy. This tendency—more obvious at 5,000 meters (fig. 47) than at 4,000 meters (fig. 
46)—was exhibited near the bathymetric feature at 30°N, 162°W, near the Aleutian Trench, 
and in the vicinity of the Hawaiian Island Rise. Throughout much of the area, circulation 
at 4,000 and 5,000 meters was opposite to that above 1,000 meters. 


Figure 48 depicts the dynamic topography adjusted to the May 1963 wind-driven 
geostrophic transport in which the Ekman divergence is neglected. Differences between 
the flow patterns shown in figure 48, which neglects Ekman divergence, and figure 47, 
which includes consideration of Ekman divergence, were not significant. 

At the northern limit of the area, near the Aleutian Islands, the dynamic topography 
suggests a boundary current in a direction opposing the flow at lower latitudes. In this 
region of high lateral friction, details of current patterns are open to some question as the 
method of computing wind-driven transport used is an approximation applicable to central 
ocean areas. 

Figures 49 through 54 depict the dynamic topography adjusted to wind-driven geo- 
strophic transport which was determined from the May 1963 mean of stresses calculated 
from twice-daily reports of pressure, using a drag coefficient of 0.0012. Results were 
similar to those shown in figures 42 through 47 (based on the higher drag coefficient and 
an alternative method of estimating mean values as discussed previously). The more 
accurate determination of wind-driven transport suggests a slightly weaker flow at the 


1,000 meter level south of 40° N. 
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3.2. Vertical Velocity 


A ; . ‘ Ow . 
The rate of change of vertical velocity with depth, zz? was determined from the ad- 
z 


justed values of meridional velocity for May, 1963, and the expression 


dw _ p 
dz ff 


ow : : 
With w assumed to be zero at the bottom, 5, Was graphically integrated between 5,000 and 
ba 


250 meters. The distributions of vertical velocity with depth are shown in figure 55. The 
maximum values determined at the four latitudes investigated varied from 3 - 10-> cm/sec 
to 14:10-° cm/sec. Though the magnitudes agree with those inferred by Knauss (1962), 
the directions do not. Knauss inferred upward velocities over the depth interval 2,500 
meters to 5,000 meters, whereas a consistent downward direction was found in this study 
at all levels below 1,000 meters. This disagreement may arise from the fact that Knauss’s 
values reflect the net circulation, while values from the present determination are asso- 
ow 
Oz 
same magnitude as the limit of sensitivity accepted for the determination. Uncertainties 
in the present determinations of vertical velocities and in the relationship of these veloci- 


ciated with a specific period. It should also be noted that the values found for — aré of the 


ties to net vertical circulation notwithstanding, it is noteworthy that the upward velocity 
found near 50°N above 1,000 meters is in accordance with flow inferred from oxygen 
profiles. In figures 5, 11, 14, and 17, the local thickening of the low-oxygen interval and 
the decrease in depth of oxygen isopleths toward the north in the depth interval 200 to 
1,000 meters suggest upward flow in this area. 


3.3 Depth of No Horizontal Motion 


In the present formulation, a component of the geostrophic velocity is zero at a depth 
where the barotropic mode is equal and opposite to the baroclinic mode. The depths of no 
zonal motion or of no meridional motion are shown as functions of latitude in figures 56 
through 62. 

If, near the depth where the baroclinic and barotropic velocities add up to zero, the 
baroclinic velocity changes slowly with depth, the determination of the depth of no motion 
will be sensitive both to observational errors in the oceanographic data and to discrepancies 
in the values for wind-driven transport. The sensitivity in the determination together with 
the paucity of data leads to uncertainty in the plots of the surfaces of no zonal or meridional 
motion. However, figures 56 through 62 do indicate that the depths of no zonal or me- 
ridional motion vary widely with time and place, and that flow at depth is frequently op- 
posite to that near the surface. 


3.4 Annual Cyele of Horizontal Velocities 


Two features of the flow at 4,000 and 5,000 meters persisted through each determina- 
tion: The well-developed westerly flow north of 45°N, and the southerly flow in the zone 
from 30° to 35°N (figs. 28, 29, 31, 32,37, 38,53, and 54). The southerly flow is in agreement 
with that inferred by Wooster and Volkmann (1960). Their data are not definitive north 
of 45°N. However, an easterly velocity can be inferred from the distributions of potential 
temperature and dissolved oxygen observed in the present study, because potential tem- 
perature increases and dissolved oxygen decreases from west to east. So it must be con- 
cluded that the velocities determined for September 1961, September 1962, and May 1963 
are not representative of the net circulation. 
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The inferred easterly net flow below 4,000 meters does not preclude the westerly 
flows determined for specific periods in this depth interval. Each of the determinations 
in the present study was based on wind-driven geostrophic transports associated with 
the weak atmospheric pressure gradients of the North Pacific summer (that is, May 
or September). Thus,.total geostrophic transport north of 45° N as found, for exam- 
ple, for the portion of Section 5 between stations 120 and 123 and 125 was small. The 
mass structure in this horizontal interval is such that a large easterly baroclinic trans- 
port existed. The velocities associated with this easterly baroclinic transport are 
strongest above 2,000 meters, decreasing to weak easterly values below 4,000 meters. 
The combination of a small geostrophic transport and large easterly baroclinic transport 
results in a large westerly barotropic transport. Above 2,000 meters, westerly barotropic 
flow is masked by the stronger easterly baroclinic flow; but below 4,000 meters barotropic 
flow is stronger than the baroclinic mode, and total flow is westerly. This result is reason- 
able for the periods associated with the present determinations. 

Because the mass structure and the associated baroclinic transport can be expected 
to be relatively constant, the annual cycle exhibited by the total wind-driven transport 
must be reflected in the barotropic mode (Fofonoff, 1961c). The baroclinic velocity has 
been shown to make but a small contribution to the flow below 4,000 meters, therefore 
a description of the annual cycle in the barotropic flow does much to describe the net 
flow below 4,000 meters. To gain a qualitative estimate of the annual cycle of the baro- 
tropic flow the transport associated with the dynamic heights through stations 120-123 
and stations 123-125 (1963) was successively adjusted to monthly values of the wind- 
driven geostrophic transport. As this facet of the investigation is qualitative at best, it 
seemed justifiable to use the most readily available tabulation of wind-driven geostrophic 
transports for the comparison. The 196] tabulations based on a drag coefficient of 0.0026 
and monthly means of the atmospheric pressure distribution were used. Figure 63 
is a plot of this determination of the annual cycle of the zonal components of velocities 
at 4,000 meters. Variations of about 1.5 cm/sec in the zonal velocities are indicated but 


the net flow is not easterly. 
The approach used for estimating the annual cycle of flow shown in figure 63 is de- 


ficient in the following respects: 

a. The monthly values used for wind-driven geostrophic transport are only approxi- 
mations of the proper values. Transport based on a drag coefficient of 0.0012 and mean 
stresses computed from twice-daily pressure reports might be significantly greater than 
those used in constructing figure 63, due to the large variance in atmospheric pressure 
which accompanies the severe, often short-period storms of the North Pacific winter. 
These larger wind-driven transports should be easterly and larger than the easterly baro- 
clinic transports, thus resulting in easterly geostrophic flow at the 4,000-meter level. 

b. The present formulation is based on the geostrophic approximation that a balance 
exists between pressure forces and Coriolis forces. It must be recognized that this 
balance is not perfectly achieved. Those time-dependent motions not satisfying the 
geostrophic approximation are not considered in the present determinations. These 
could make a significant contribution to the net flow which has been inferred from the 
distribution of potential temperature and dissolved oxygen. 

c. The present investigation is based on the assumption that the total wind-driven 
geostrophic transport and the barotropic transport respond significantly only to disturb- 
ances of the one-month period considered. Veronis and Stommel (1956) indicated that 
the barotropic response might be significant for disturbances with periods of from one 


week to two months. 
d. In the comparisons supporting figure 63, it is assumed that baroclinic transport is 


constant and can be accurately described by dynamic height computations based on 
oceanographic data from May 1963. Obviously, there may be a significant variation in 
the mass structure and thus in the baroclinic transport through the section. 
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e. Although the baroclinic mode has been shown to be small in comparison to the 
barotropic mode below 4,000 meters during three specific periods, it cannot be assumed 
that the annual mean baroclinic flow is small in comparison to the annual mean barotropic 
flow. Veronis and Stommel (1956) indicated that net circulation is associated with the 
baroclinic mode. It can be anticipated that an accurate description of the annual mean 
baroclinic transport through stations 120-123 and 123-125 might indicate a net easterly 
flow at the 4,000 meter level. 
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4. Summary and Conclusions 


A method has been derived for matching the geostrophic transport, calculated from 
dynamic heights, with wind-driven transport for the same period. Geostrophic velocities 
at 5,000 meters and above were determined covering a portion of the north central Pacific 
Ocean for the months of September 1961, September 1962, and May 1963. These ve- 
locities are reported as plan views of dynamic heights. The determined velocities were 
generally zonal for the periods of observation. At the 5,000-meter level values on the 
order of 1 cm/sec were found. Over most of the area at these times, the sense of the 
velocity changed within the water column, so that the current direction at the 5,000-meter 
level was opposite to that of the near-surface levels. 

The effects of varying the drag coefficient used in determining wind stress and of 
alternate methods of determining the mean stress distribution for the one-month period 
of consideration was investigated. The use of a drag coefficient of 0.0012 and values for 
stresses determined as the mean of stresses calculated from twice-daily reports of pressure 
distribution seem to give more realistic transports than those obtained from other com- 
binations of drag coefficients and pressures. 

The horizontal distribution of potential temperature and dissolved oxygen indicate 
that the net flow at 4,000 meters at about 48°N must be easterly. Such a flow is ina direc- 
tion contrary to the flow found at this location by the transport method. This has been 
accounted for by examining the relative importance of the baroclinic and barotropic modes 
at this level. The flow at 4,000 meters is mainly barotropic, and as the barotropic flow 
must be periodic (Veronis and Stommel, 1956) it is concluded that the calculated deep 
flow is a periodic barotropic flow. The calculated magnitude of deep velocities (about 
1 cm/sec) is in good agreement with Veronis and Stommel’s theoretical value. 

Vertical velocities are determined by graphically integrating o determined as a func- 
tion of meridional velocity. Values as high as 10-4 cm/sec were calculated, but these 
are at the level of significance of observational errors. Thus, the results shown may re- 
flect the range of variation in vertical velocity, but details of the variation with depth cannot 
be validated. 

For each period for which horizontal velocities were computed, plots were made of 
the surfaces of no zonal or meridional motion. The determination of the depth of no mo- 
tion can be very sensitive to observational errors and inaccuracies in the computed wind- 
driven transport; consequently, the plots are subject to considerable error but do indicate 
that the depth of no zonal or meridional motion varies widely with place and time. 

The assumption that terms due to the interaction of flow with the bottom can be ig- 
nored was examined. The data suggest that bottom flow parallels isobaths and therefore, 
the bottom and slope terms in the integrated momentum equations become very small. 
However, determinations in greater detail over an area in which the bathymetry is well 
defined would be needed to show conclusively that these terms are negligible. 

Because the velocities determined near the bottom appear to be periodic, comparison 
with net velocities as inferred from the horizontal distribution of variables (Wooster and 
Volkmann, 1960) or by carbon-14 dating or heat flow considerations (Knauss, 1962) is in- 
appropriate. Perhaps the best corroboration of the deep flow determined in the present 
study is in the tendency of the water to flow parallel to isobaths. 

It is concluded that the method employed here is a useful tool for the study of circula- 
tion in central ocean areas. To best realize the possibilities of the method, oceanographic 
data should be collected so that no additional limitations are placed on the results. A 
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grid of stations forming sections along meridians and parallels approximately 300 kilo- 
meters in length would allow full realization of the capabilities of the method. In addition 
to a more detailed determination of the flow patterns, such a grid would provide means 
better to determine the validity of deep observations and to assess the importance of in- 
ternal structure to the deep flow. Observations in the grid should extend to within a 
few hundred meters of the bottom. 
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Appendix I. 


Physical Expressions for the Problem 


Sverdrup (1947), Munk (1950), and Fofonoff (1962) have developed expressions relat- 
ing total transport to the curl of the wind stress. All three formulations were predicated 
on a velocity structure such that velocity vanished well above the bottom. Thus, it was 
not necessary to consider the effects of bottom topography on the flow. In the present 
problem, portions of these existing solutions applicable to central regions are employed. 
Approximations and assumptions are made pertaining to the effects of bottom topography 
on the solution for transport. 

The development is based on Fofonoff’s (1962) treatment of Munk’s (1950) wind- 
driven flow formulation which is in turn, developed from the momentum equations for 
steady motion. In a coordinate system where x is positive eastward, y is positive north- 
ward, and z is positive upward and assuming incompressibility and negligible molecular 
stress terms, the following equations are obtained: 


dpu*  dpuv , dpuw oP ‘5 Ou 
Oe OO pfu =—=— + pV? aS 
SM TL enor Utes 
dpuv , dpv? _ dpvw oP z d*v 
ae a ee + pfu=——+F pV? 
FEMA MLNS IM camer IN ee 
BY la, b 
ae (1a, b, c) 
and continuity is expressed by 
d(pu) , d(pv) , (pw) __g 2) 


Ox oy Oz 


Assuming terms with vertical velocities are negligible, the vertically integrated equa- 
tions of motion can be written as 


7 dpu? ” dpuv i = 2 aP i (3+ u ane: = a 
I ae dz+ 4 ([ dz—f | pvdz= ae dz+ An fl DA eer dz+ 

7 dpuv dpv" = aP i eS v an =) I Ov 
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where f is the Coriolis parameter, p is mean density, and Ay and Ay are the horizontal and 
vertical kinematic eddy viscosities, respectively. 

Convenient terms for the integrals are derived by the following definitions and argu- 
ments. First, define the components of horizontal momentum transport 


w= |" puvdz, =|" pu*dz: (3a, b) 
a 2p 


and the components of mass transport, U and V, 


n n 
u={ pudz, v={ pvdz, 
2p 2p 


2 GS * Die 12 Ph a (4a, b) 
xX x 


and, since Py and 07/dx are both very small, 


pea-2 [7 Paes Poe, Pea-Z[ Pdz+ P= 
2b di) 


Next, define the potential energy of a column 


n 
Ep =| Pdz 
*p 


so that when equation (6) is substituted in (5a, b) 


n 
{ La eels ag PE tes, 
zp) OX Ox Ox zy OY oy ay 


From equation (1c) 
éz=—— aP 
& 


and on substitution into (6) 


(5a, b) 


(6) 


(7a, b) 


(8) 


Letting the specific volume, a=ao+6 wherein q is the specific volume of water of sa- 


linity 35°/o0, temperature 0° C, and pressure P, 


P P 
=| " PaglP + | ” PsdP 
& Jo 2.0) 


Fofonoff (1962) defines the last term in (9) as the anomaly of potential energy, 


-1{” PoaP. 
& Jo 


Using the definition in the first term on the left of equation (7a), 


ceo ONE ae I, ’ PaudP |, 


ox oO pes 


and by application of the chain rule 


dEy_ ax , Py a 


1 
dx ax ax raat  PouaP | 


(10) 


Since the bracketed portion of the last term is a function of the bottom pressure only, 


OE _ OX AP» Pra Ps) 
Ox Ox Ox g& 
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This, and the similar expression for any substituted in equations (7a, b) give 
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Finally, it can be shown that the final terms 


a O-a 
Ay dz=Tsx — Tox, 
z 


2 
> Oz 
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[. Ay a dz= Tsy — Toy, (12a, b) 


where 7s, represents the stress per unit surface area in the x direction and 7), is a similar 
term for stress at the bottom. By substituting equations (3a, b), (4a, b), (lla, b), and 
(12a, b) into (1a, b), 


uz duv dx AP» Pyao(Py) 02 : 
oon =——4 h PY ZU) Sc 
ae ay iV aE ax Z Py ae AnV U = Uti = Uta 
duv , av? Ox APs Proo(P, a aoe 
sat 5 +fU= ie a Ze ») _p, ay AVN + Tey — Toy. (13a, b) 


The bottom pressure and bottom slope terms in (13a, b) can be further simplified as follows: 


= OP, dP: |e OZ 
[Bix Pap ie dz, =—+ e : 
z i «| : OX eh, Ox dz— go Ox 


When z=z the second term on the right is eliminated so that 


OP» Ns eS) 5 O20. 
dx NOxyeae ax 


With this substitution the terms from (13a) become 


Pyao(Po) () 02 OZ 
ae sacrament || Sa et + Py iLO aa Rare 
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It is assumed that ao(P,)=a, and the bottom pressure and slope terms of (13) may be 
represented by the single term EHO (=) 
& dx 2=25 

To be useful for the present purposes it is necessary to neglect additional terms in 
(13a, b). Assume (a) that field acceleration terms are negligible; (b) that the velocities 
near the bottom are very small so that 7 is much smaller than 7;, so the vertical eddy vis- 
cosity terms when integrated over the column are closely approximated by 7;; and (c) that 
in a central region of the ocean horizontal friction terms are negligible (Munk, 1950). With 
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these assumptions, the transport equations may be written: 


P, 
fv =—2X_ om i + Tsx, 
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There are no existing methods of determining open ocean transports with sufficient pre- 
cision to account for mass exchange at the sea surface by continuity methods. To a 
sufficient degree of accuracy the divergence of the transport is expressed by 


aU , av _ 


ax day o (15) 


Equations (14a, b) and (15) are the basic ones for this problem. 


After Fofonoff (1962), the components U and V of the transport can be separated into 
barotropic, baroclinic, and Ekman modes: 


U=U,+ U,+ Ur 
(16a, b) 
V=Vi+Va+Ve 


The barotropic mode is interpreted as a velocity, vp, constant throughout the water column, 
so that 


Pe _ (aP, 
fur aa ( Oy )_, fry ak ls I 
n n 
i I bibs. VE I Po (17a, b) 
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since uw and vp are not functions of z 


Ur=w |” pds Vo=m | pd. 
2b z 


By the hydrostatic equation 


P P 
= | "dP, =Vy=2 | ° aP; 
& Jp, 
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U,= uoPo V,= wPo 
& & 
By substituting for u» and vp in equations (17a, b) 
_ fy, = o0Po 2) — Po GS 
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The baroclinic mode is interpreted as that part of the transport that is a consequence of 
pressure gradients arising from field variations in the density structure of the water. The 
velocities associated with this mode of transport are the relative velocities resulting from 
dynamic height computations. This mode of transport can be related to the pressure 
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terms on the right side of equations (14a, b) as follows: 


oP 
Pf (vy + vg) = Foe 
and integrating 
oP 
(sie || Se 
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Substituting for VY; from (18b), 
aP, (AP: 3 {vole 
& ( Ox ema Ox ce: 


Then using (11a) and the arguments following (13a, b) regarding bottom pressure terms, 


oP (Pz _ ay , aPs (AP: 
& feel ee ae & eal 7 (19) 


Identical pressure terms can be eliminated from the two sides of this expression. This, 
together with a similar development for U, yields 


0 Ox 
U, ee V, =: 
LU, ay Wim (20a, b) 


The components of the Ekman transport are defined by the equations 


LUE = Tey, Vn=—Tsz- (2la, b) 


Although the sum of the divergences of the three modes must be zero (15), the divergences 
of each mode taken separately need not be zero. Expressions for the individual diver- 
gences can be derived by cross differentiating the equations that define the barotropic, 
baroclinic, and Ekman transports—(18a, b), (20a, b), and (21a, b). Cross differentiating 
leads to 


au, OV, = O02 02 
A at ai ) + BV, = pf (ws + vp a (22) 
Aa a1) + AV, =(), (23) 
0 Oy 
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In equation (22) p’ is closely approximated by the mean density of the column. 

The right side of equation (22) can be neglected when the gradient of the bottom is 
sufficiently small or when flow along the bottom is nearly normal to the slope of the bottom. 
(Applicability of this approximation is examined in section 2.2.) With this approximation, 
on adding equations (22), (23), and (24) 
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and on substituting for transports from equation (16a, b), and noting the continuity condi- 
tions expressed by (15) 


OTsy OTs 
1 ee ee ‘ 
B ax ay (25) 
This is the expression derived by Sverdrup (1947). 
Equation (15) allows the introduction of a transport function, ws, such that 
aw aw 
U=—; V=——- y 
ay ae (26a, b) 
Substituting from (26b) into (25) and integrating yields 
=| (ee aa diene. (27) 
B 
Using conditions at the eastern boundary, the constant 
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ar le 
on combining with (27) 
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This solution (Fofonoff, 1962) does not provide for an eastern boundary layer. Fofonoff 
has justified the elimination of an eastern boundary layer by writing Munk’s (1950) general 
solution for mass transport in a nondimensional form and simplifying the equation on the 
basis of relative orders of magnitude. For the area under consideration (the central 
North Pacific Ocean) provision for an eastern boundary layer is not necessary. 

As mentioned above Sverdrup (1947), Munk (1950), and Fofonoff (1962) all specified 
that the transport was made up of baroclinic and Ekman modes only; the problem of 
treating convergence due to bathymetric effects on the barotropic flow thus was not 
pertinent. 

Fofonoff (1960) used equation (28) (in geographic coordinates) as the basis for his 
mass transport computations in the Pacific. Using monthly means of sea-level atmos- 
pheric pressure (from the records of the U.S. Weather Bureau) he successively computed: 

geostrophic wind velocities, 

surface wind velocities, 

wind stress on the surface, 

meridional and zonal Ekman transport, Vz and Uz, 
meridional total transport, V’, 

integrated total transport, W, and 

integrated geostrophic transport, wy. 


The pressure distribution and each of the components or modes of transport are tabulated 
on alternate five-degree grid points for the Pacific Ocean north of latitude 20°N. Tabu- 
lations have been published for each month during the period January 1950 through May 
1962 (Fofonoff 1960, 1961: Fofonoff and Ross, 1962: Fofonoff and Dobson, 1963). Similar 
tabulations have been made for months of interest during the period June 1962 through 
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June 1963 by the Applied Mathematics staff of the Department of Oceanography, Univer- 
sity of Washington. 
Transformed into geographic coordinates the quantities of interest are: 


Ve=-—2, Up =*8, 
ar a chant (29a, b) 
__ (dT  O(T. cos 2) I 
V : 
ee od BR cos '% (30) 
The transport function, w, is defined by 
1 ow 1 ow 
——— Vy =—-—— —. 
R dd R cos dX (31a, b) 
Equation (31b) is integrated from the eastern shoreline, 
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Fofonoff (1960) tabulates the geostrophic transport as a transport function, j. He notes 
this this procedure is not entirely correct, since the geostrophic transport is not divergence- 
less. He obtains his values for W from the expression 


w=. [= (Be Te ‘| dn, (33) 


that is, the difference V—V, integrated over longitude. Following this procedure the 
meridional geostrophic transport is accurately represented by J; however, the zonal 
geostrophic transport is incorrectly represented since it reflects the effects of neglecting 
divergence of the Ekman transport. 

Both the transport function for total transport and that for geostrophic transport can 
depict zonal transport in error as a consequence of using the grid point nearest the eastern 
boundary as the origin for integration. However, in the present use of Fofonoff’s tabu- 
lated transports as reference values, results are insensitive to the wind-driven transport 
to such a degree that errors caused by unaccounted-for zonal transport near the eastern 
boundary will not greatly influence the final values. 

With the values tabulated by Fofonoff (1960), baroclinic transport computed from 
oceanographic station data, and equation (16a, b), it is possible to obtain barotropic trans- 
port through a given section: 


U,»=U—Ur—U,y, 
(16a, b) 
V,=V— Vira 


In the present problem the intervals of interest are determined by the locations of the 
oceanographic stations. The deep stations used are located so as to form either zonal 
or meridional sections. In equations (16a, b) the first terms on the right can be extracted 
from Fofonoff’s (1960) tabulated transport functions for total transport by finite differ- 
ences; the second terms, Ekman transports, can be integrated over the interval from the 
tabulated values; and the third term can be taken from the oceanographic data. Equations 
(17a, b) can be closely approximated by 


Up = PmUph 5 
(34) 
Vy = Pmvyh. 
25 


798-230 O-66—3 


And so it is possible to solve for the barotropic velocities up, vp 
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ey ae Se eee 
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These barotropic velocities, constant throughout the column, can be combined with baro- 
clinic velocities based on observations over the entire depth of the column to obtain 
absolute geostrophic velocities. 
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Appendix II. 


A Numerical Example 


Based on equations (35a, b), velocity profiles are determined for the meridional and 
zonal sections defined by oceanographic stations 29, 32, 59, and 63, September 1961. 

The first step in the determination is to extract the geostrophic mode of transport 
for the section of interest. This transport is included in equations (35a, b) as the first 
and last terms in the numerator on the right side. Values of the total transport function 
along the meridians of interest were plotted against latitude (fig. 64). The stations are 
indicated on the plots of W, and the total transport through the sections between the indi- 
cated stations is taken by difference: 


Section 

63-59 (Ws, — ¥s,) =[ — 187 —(— 112)] - 10° =— 250 - 10* metric tons/sec 
29-32 (ena) | 92 ( a3) |i 10® —— 790) 102 metricitons/sec 
63-29 QU) lls (12) 05 —— 105-10" metric tons/sec 
59-32 (Wh, — vx) =[- 192 —(— 137)] - 10°=— 550 - 10* metric tons/sec 


{t can be seen that continuity requirements for the total flow are satisfied by the results 
of this procedure. Plots of zonal and meridional Ekman transport along meridians and 
parallels are shown as figures 65, 66, and 67. Again the stations are indicated on the plots: 
summations to obtain Ekman transport through the sections are performed using Simpson’s 
rule, or the trapezoid rule: 


Section 
63-59 SU; Rap— | a0 

= 3.2 - 104 metric tons/sec 
Zone 3 Rap=| DUES 9S) SHS 8) 


=4.8 - 104 metric tons/sec 


63-29 SVR cos {q- ee eee | 10 


3 3 
= 64.1 - 104 metric tons/sec 
1520(2 + (—2. 
59-32 YVsR cos ddd= | - 10 
=—(0.2 - 10* metric tons/sec 


The wind-driven geostrophic transport for a given section is the difference between the 
total and Ekman transports: 


Section 

63-59 = (— 250—3.2) . 10! =— 253 - 104 metric tons/sec 
29-32 = (— 790— 4.8) - 104 =~—795 - 104 metric tons/sec 
63-29 = (—10—64.1) - 10!=— 74 - 104 metric tons/sec 


59-32. =(—550—(—0.2)) - 10* ~— 550 - 104 metric tons/sec 


Because the distribution of horizontal velocities is not highly sensitive to this geostropnic 
transport, values for it can be extracted directly from tabulations of the geostrophic 
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transport stream function, j. (As noted in appendix I, this procedure eliminates from 
consideration the Ekman divergence.) Plots of the geostrophic transport function similar 
to those for the total transport function are made (fig. 68) and the geostrophic transport is 
extracted by difference. A comparison of alternate values of geostrophic transport for 
the four sections considered is given below. 


Section Transport considering Transport from geostrophic 
Ekman divergence transport function 
63-59 — 253-10! metric tons/sec — 20-10! metric tons/sec 
29-32 — 795-104 metric tons/sec — 310-10! metric tons/sec 
63-29 — 74-104 metric tons/sec — 250-10! metric tons/sec 
59=32 — 550-104 metric tons/sec — 540-10? metric tons/sec 


The dynamic topographies arising from a comparison of transport from dynamic heights 
adjusted to each of these sets of wind-driven transport are shown on figures 23 through 29. 
For the balance of the numerical example the transport including the effects of diver- 
gence is used. 

The final terms in equations (35a, b) express mass transport through the sections due 
to internal density structure. To obtain this baroclinic transport, data from each station 
is submitted to a machine interpolation program which provides the following values at 
standard depths: 

temperature, 

salinity, 

sigma-t, 

specific volume anomaly, 6, 

geopotential anomaly, AD, referred to the surface, 
anomaly of potential energy, x, 

and values for other variables represented in the field data. 

For convenience in the computations advantage is taken of the near numerical equality 
between mass transport and volume transports, 7 and 7;. The mass transports deter- 
mined from the curl of the wind stress are entered in subsequent comparisons as though 
they were volume transports. The final terms in (35a, b) are determined from differences 
in geopotential anomalies. Equations (35a, b), reflecting these near-equalities, become: 


AD, — ADa 
p—7.—(SBo 8h) 
hi — f 
2 hRAd 
rT. (2 z a) 
i AR cos @AX (36a, b) 


Next, values are determined for the length-of-section terms, R cos @AA and RAd, 


1 . ore 
and for —, the reciprocal of the Coriolis parameter: 


if 
Section Mid-latitude Length (km) jlo sec) 
63-59 34°50’ 480 0.1201 
29-32 34°33’ 485 0.1209 
63-29 O2Role 1613 0.1276 
H9=52 30102) 1520 0.1141 
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The depth of the column, h, was approximated as 5,000 meters, and geopotential 
anomalies from the station data were referred to the same depth. Observations ex- 
tended to at least 5,000 meters at most of the stations used in the program; at stations where 
they did not, values of temperature and salinity were selected to extend the profiles to 
5,000 meters. These values were selected by extrapolating profiles at the station con- 
sidered and by interpolating between appropriate values at adjacent stations. In the 
selection care was taken to avoid creating artificial structure in the water mass. 

Values of the wind-driven transport, geopotential anomalies referred to 5,000 meters, 
the distance between stations, and the reciprocal of the Coriolis parameter are submitted 
to a second program. This program calculates velocities and geopotential anomalies 
at each standard depth. These velocities and anomalies reflect the variation with depth 
associated with the structure of the water column and are adjusted to describe a trans- 
port through the section equal to the geostrophic mode of transport from the curl of wind 
stress — that is, equations (36a, b) are satisfied. Tables 2 through 5 show the data printouts 
from this program for the four sections considered. An explanation of the quantities 
shown in the headings and columns of each table follows. 


Quantities in the headings are: 

STATIONS, identifying numbers of oceanographic stations 

L, length of section in kilometers 

COR, reciprocal of the Coriolis parameter, multiplied by 10-° 

TRANSPORT, the wind-driven geostrophic transport in 104 m3/second 

N, the number of standard depths considered 

SUMVZ, transport in 10! m®/second as obtained from differences of dynamic heights, 
referred to 5,000 meters; the sum of the values in column VVZ below 

CORR, a correction factor in volume transport per meter of depth in 10” m?/second, 
obtained by: 

TRANSPORT — SUMVZ 


CORR= h , that is 


AD, — ADa 


1 fig 7 


CORR= 


SUMVZ*, acheck of computations. The value is identical with that for TRANSPORT. 
Quantities in the columns are: 

DEPTH, standard depths 

A, geopotential anomalies in dynamic meters for station A, referred to 5,000 meters; 
station A is the southern station of meridional sections and the eastern station of 
the pairs along parallels of latitude. 

B, geopotential anomalies for station B 

C, the difference, column A minus column B at each standard depth 

VZ, volume transport per meter of depth in 10? m?/second applicable to each depth. 
This function of volume transport is obtained by multiplying the difference in geo- 
potential anomalies (column C) by the reciprocal of the Coriolis parameter (COR in 
heading). 

VVZ, baroclinic transport in 10! m?/second for each depth interval, obtained by 
numerical integration of VZ with respect to depth 

VCZ, an adjusted volume per second per meter of depth applicable to each depth. 
VCZ is obtained by adding VZ at each depth to the correction, CORR, as given in 
the heading of the table. 

V*, the velocity at each standard depth corresponding to the adjusted transport values, 
VCZ. These velocities, obtained in cm/second, are the absolute geostrophic 
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velocities. The value of V* at the 5,000-meter level corresponds to up or vp in 
equations (36a, b). 

C*, adjusted differences in geopotential anomaly for the section. These adjusted 
differences are accepted as reflections of the absolute geopotential 


topography. 
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FIGURE 2.—Stations used in determination of transport. 
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FIGURE 3.— Temperature (°C) Section 1, September 1961. 
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FicurE 4.—Salinity (Yoo) Section 1, September 1961. 
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FicurE 5.—Dissolved oxygen (ml/L) Section 1, September 1961. 
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FIGURE 6.— Temperature (°C) Section 2, September 1961. 
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FicureE 7.—Salinity (00) Section 2, September 1961. 
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FIGURE 11.—Dissolved oxygen (ml/L) Section 3, September 1962. 
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FIGURE 12.—Temperature (°C) Section 4, September 1963. 
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FIGURE 13.—Salinity (oo) Section 4, September 1963. 
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FIGURE 23.—Topography in dynamic meters, 250-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence con- Oe 
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FIGURE 24.—Topography in dynamic meters, 500-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, September 1961. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 
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FicuRE 25.—Topography in dynamic meters, 1,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, September 1961. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 
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FIGURE 26.— Topography in dynamic meters, 2,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, September 1961. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 
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FIGURE 27.— Topography in dynamic meters, 4,000-m level adjusted 
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FIGURE 28.— Topography in dymanic meters, 5,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, September 1961. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 
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FicuRE 29.—Topography in dynamic meters, 5,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence not con- 
sidered, September 1961. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 
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FIGURE 30.—Topography in dynamic meters, 2,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence not con- 
sidered, August 1961. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 
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FIGURE 31.—Topography in dynamic meters, 4,000-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence not con- e oe 
sidered, August 1961. Drag coefficient 0.0026; stresses from g 
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FIGURE 32.—Topography in dynamic meters, 5,000-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence not con- re 0.8 
sidered, August 1961. Drag coefficient 0.0026; stresses from ~ 
mean monthly pressures. 2 es 
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FIGURE 33.—Topography in dynamic meters, 250-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, September 1962. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 


69 


VELOCITY, CM /SEC 


20° 30° 40° 50° 
VELOCITY SCALE FOR .O! METER INTERVAL 


FIGURE 34.—Topography in dynamic meters, 500-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, September 1962. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 
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FIGURE 35.—Topography in dynamic meters, 1,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, September 1962. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 
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FIGURE 36.— Topography in dynamic meters, 2,000-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence con- ey Sie 
sidered, September 1962. Drag coefficient 0.0026; stresses from we 
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FIGURE 37.— Topography in dynamic meters, 4,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, September 1962. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 
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FIGURE 38.— Topography in dynamic meters, 5,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, September 1962. Drag coefficient 0.0026; stresses from 
mean monthly pressures. 
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FIGURE 39.—Topography in dynamic meters, 5,000-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence not con- ° OS 
sidered, September 1962. Drag coefficient 0.0026; meaned Q 
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FIGURE 41.— Topography in dynamic meters, 5,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence not 
considered, September 1962. Drag coefficient 0.0012; meaned 
stresses from twice-daily reports of pressure. 
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FIGURE 42.—Topography in dynamic meters, 250-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence con- i oe 
sidered, May 1963. Drag coefficient 0.0026; stresses from mean xz 
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FicuRE 43.—Topography in dynamic meters, 500-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence con- = eke 
sidered, May 1963. Drag coefficient 0.0026; stresses from mean g 
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FicureE 44.—Topography in dynamic meters, 1,000-m level adjusted 
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FIGURE 45.— Topography in dynamic meters, 2,000-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence con- o O89 
sidered, May 1963. Drag coefficient 0.0026; stresses from mean x 
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FIGURE 46.—Topography in dynamic meters, 4,000-m level adjusted 


e F : fi 0.5 
to wind-driven geostrophic transport, Ekman divergence con- 8 
sidered, May 1963. Drag coefficient 0.0026; stresses from mean ~ 
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FIGURE 47.—Topography in dynamic meters, 5,000-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence con- ° 0.5 
sidered, May 1963. Drag coefficient 0.0026; stresses from mean 2 
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FIGURE 48.—Topography in dynamic meters, 5,000-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence not con- oro: 
sidered, May 1963. Drag coefficient 0.0026; stresses from mean g 
monthly pressures. 3 10 
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FicurE 49.—Topography in dynamic meters, 250-m level adjusted 
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to wind-driven geostrophic transport, Ekman divergence con- 3 
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FIGURE 50.—Topography in dynamic meters, 500-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 


sidered, May 1963. 
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FIGURE 51.—Topography in dynamic meters, 1,000-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence con- i o:8 
sidered, May 1963. Drag coefficient 0.0012; meaned stresses from g 
twice-daily reports of pressure. 5 16 
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FIGURE 52.—Topography in dynamic meters, 2,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, May 1963. Drag coefficient 0.0012; meaned stresses from 
twice-daily reports of pressure. 
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FIGURE 53.— Topography in dynamic meters, 4,000-m level adjusted 


to wind-driven geostrophic transport, Ekman divergence con- ° Of) 
sidered, May 1963. Drag coefficient 0.0012; meaned stresses from R 
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FIGURE 54.— Topography in dynamic meters, 5,000-m level adjusted 
to wind-driven geostrophic transport, Ekman divergence con- 
sidered, May 1963. Drag coefficient 0.0012; meaned stresses from 
twice-daily reports of pressure. 
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